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Abstract 

We study the asymptotic behaviour of resistance scaling and fluctuation of resistance that 
give rise to flicker noise in an n-simplex lattice. We propose a simple method to calculate 
the resistance scaling and give a closed-form formula to calculate the exponent, {3l, associated 
with resistance scaling, for any n. Using current cumulant method we calculate the exact noise 
exponent for ra-simplex lattices. 

1 Introduction 

In recent years considerable attention has been devoted to studying the properties of disordered 
systems with the hope of understanding percolative phenomena. Key to several such approaches 
has been the concept of randomness and also of frustration [1-12]. However, many of the patterns 
we encounter in nature are not random but self-similar and scale invariant [13-14]. For instance, the 
complicated and scale invariant structures that occur when a solid mixture evolves via an aggregation 
process [15]. To understand such systems the concept of fractals has been found to be very useful. 
Fractals are scale invariant objects that may be considered as intermediate lattices between regular 
and random (disordered) lattices [13-14,16-17]. Such a fractal lattice describes a class of random 
systems where the consequence of the loss of translational invariance of a lattice can be studied in 
detail. Additionally, resulting from their dialational symmetry, statistical, mechanical and transport 
problems are solvable; hence the attraction of the model in such studies [14] . 

In this paper we consider a particular class of fractal known as the n-simplex lattices, to model 
various properties of inhomogeneous materials [16-17]. The lattice is defined recursively. The map 
of the zero-order truncated n-simplex lattice is a complete set of (n+1) points. The map of the 
(r+l)th order n-simplex lattice is obtained by replacing each of the lattice points of the rth order 
map by the entire rth order map. Each of the resulting n points is connected to one of the lines 
connecting the original rth order vertices. The fractal and spectral dimensions of this lattice are 
given by: 

and 

d -ln(n + 2)- {Z) 
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The lattices with n > 3 are of particular interest as they provide a family of fractals in which d 
varies with n leaving d almost constant. 

In order to understand how resistance scales with the size of the system, in a homogeneous 
system, we study the distribution of currents in a network modeled by an n-simplex lattice. We 
consider each bond, where bond refers to a line joining two lattice points, of the zero-order network 
as a unit resistor offering resistance R. A unit current enters the network at one of the external 
nodes and leaves through another, the rest of nodes being left open. It is known that the distribution 
of currents in such a network is found to be multifractal [18] in the sense that different moments of 
the distribution scale with different exponents. 

For resistance scaling analysis two methods may be adopted, cither (a) to obtain the distribution 
of current over the entire network and measure the energy dissipated in the system or (b) to simplify 
the network and obtain a closed-form solution. This second method has been used rigorously by 
us for resistance scaling and the results obtained by this method match those from the current 
distribution method. The moments of the current in a n-simplex are: 

s a r {h,h, /„) = £|/ P r, (3) 

p 

where I p is the current in the pth bond and p goes from 1 to n(n — l)/2; S® is the cumulant for an 
arbitrary exponent a. The currents flowing in at the external nodes of a n-simplex are represented 

by h, h, I n respectively (see figure 1) with the condition h + h + I n = 0. A scaling 

factor independent of I\, I2 I n can then be defined as: 

w x S?+i{h,h, In) ,.s 

b?(h,l2, In) 

Note that A (a) is related to the fractal scaling exponent D(a). For a fractal with a resis- 
tance scaling parameter of 2, the rth generation length scales as L r = 2 r . Using the definitions 
S?(h,h, In) = L? {a) and S«{h,h, /„) oc X r (a), we get: 

D(a) = (5) 

The case a = determines the fractal dimension of the simplex because A(0) is simply the ratio 
of number of bonds in successive order of the n-simplex lattice; a = 2 measures the heat loss in the 
network and gives resistance scaling. It has been shown [19-20] that: 

R{L)~L~fo (L»l) (6) 

where /?; is an exponent controlling the transport properties. 

In disordered material, the elastic scattering of the carriers at impurities leads to the random 
conductance or resistance fluctuation. The fluctuation arises from the interference of the scattered 
waves, and they are random. The magnitude of the resistance noise spectrum (flicker noise 1/f) 
depends on a new exponent, b, pertaining to the fractal lattice. This exponent (corresponding to 
a = 4) is a member of infinite number of exponents required to characterize the fractal lattice [18]. 
The exact reason as to why this fluctuation occurs is unknown though it is believed that it appears in 
response to changes in many extrinsic parameters such as the carrier density, the applied measuring 
current, external electric fields and external magnetic fields. The spectrum of resistance fluctuation 
is given by 

S R (w) = J e iwt < R(t)R(0) > dt. (7) 
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The exponent b associated with the scaling behaviour of normalized noise is given by [19]: 

PR = U/R 2 ~ L-» (L » 1) (8) 



As long as each bond resistance fluctuates independently with the same spectrum, the explicit 
frequency dependence can be discarded. The upper and lower bounds of b [19-20] are given by: 

/3 L <b<d. (9) 

The paper is organised as follows: In Section 2 we derive a closed-form solution to calculate (3l 
for any value of n. In Section 3 we use the current cumulant method to calculate the noise exponent 
for n-simplex. The paper ends with a brief discussion on the bounds proposed and comparison with 
our results with experimental data. 



2 Calculation of (3l associated with resistance scaling 

In this section we propose a simple method of calculating /3j for any n-simplex. Consider a fixed 
current I\ entering at one of the external nodes of the lattice, and leaving from another, all the 
remaining external nodes being left open (I 2 = I3 = ...I n = 0). We calculate the equivalent resistance 
between these two external nodes and establish a recursion relation between rth and (r+l)th order 
lattices and use the Real Space Renormalization Group Technique to find the exponents [18,21-22]. 

^From the symmetry properties of the simplex, it is apparent that all (n-2) external nodes apart 
from those through which current enters and leaves are at equipotential. Redrawing just those bonds 
through which currents flow, we have (n-1) parallel paths for current to flow. Of these paths, one 
offers unit resistance and each of the others offer twice the unit resistance (since they include two 
resistances in series). Hence the equivalent resistance is given by: 

lk = R- + [ jR- + ^R- + --- {n - 2)tCim8] (10) 

where R is the unit resistance and the square bracket contains exactly (n-2) identical terms. This 
directly leads to: 

Re = — . (11) 

n 

Now if we consider a star of ^-branches, each offering a resistance of R/n, the effective resistance 
between any two external nodes through which current flows will be 2R/n as they are in series and 
all other nodes being left open. It is then straight forward to show using these transformation for 
n-simplex lattice that the following scaling holds good: 

,/From equation (11) we know that for any n-simplex the equivalent resistance of first order is: 

2R , . 

Rei = — , 13 

n 

combined with equation (12) gives the equivalent resistance of rth order as: 

^ (14) 
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Thus we see how, by merely knowing the simplex one can calculate the equivalent resistance of 
any iteration. No long winded applications of Kirchoff 's Laws are required to obtain the resistance 
scaling. The exponent (3l is related to A(2) by: 

ft-*™. <»> 

3 Calculation of the Flicker noise exponent on n-simplex 

It has been shown that the A th moment of current distribution is associated with the noise exponent 
[18]. 

Assuming that the cumulant S^{I\,l2, I n ) can be expressed as a homogeneous polynomial of 

degree 4, the most general polynomial is a linear combination of P-j 4 , P 2 P 2 , P1P3, P4 and P|. These 
polynomials are defined as 

Pi = h+h + h + I n 

I'.' = /;; 

= / 1 3 + / 2 3 + / 3 3 + /;: and 

pa = I n 



However in present case P\ = due to current conservation. Hence Sf{I\,l2, I n ) can be 

written as 

S${h,h, In) = A r Pi{I U h, /„) + B r Pi{I U h, I n ) (16) 

The next step is to determine S^_ 1 (Ii,I 2 , In)- To establish a recursion relation between 

Sf(Ii,l2, I n ) and S'*_ 1 (/i,72, In) we obtained current distribution in an n-simplex lattice at 

each node. It is easy to see that the current distribution at each node is 



E 

fe=i 



h + 



n 



by current conservation. In figure 1 we have shown the current along each bond. Therefore, we can 
write 

S^Iuh, /„) = S r 4 -iOO + + + S^in) 

where S r -\(I), S r -i(H), are the current cumulants of (r — l)th order of n-simplex lattice. /, 

II, represents shaded region in figure 1. Above equation can be expressed as 



Sr(h,h, In) — ^r-lPUhjI^, In) + B r _ 1 P i (I 1 , I 2 , I n ) 

+ A r -iPi(Ii,l2, I n ) + B r _iP2(h,l2, In) 



(17) 



which establish the transformation relation between r and (r — l)th order. Comparing equations 
(16) and (17) we obtain the recursion relation between the n-simplex lattices of the r and (r — 1) th 
order. 

' A r \ 1 /(n 3 + 2)n n 2 (n + if \ f A r 
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(2n + 3)n J \B. 



i r _i 
? r _i 



(18) 
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The eigenvalues corresponding to the transformation matrix for the n-simplex lattice are given by 



, _ _ (n 3 + 2n + 5)n ± n(n + 1 )Vn 4 - 2n 3 - n 2 + 2n + 25 
K [a -A)- — 

and the fractal scaling exponent corresponding to largest eigenvalue is 

£(q = 4) = lnA " (a = 4) . 
v ; In 2 

4 Discussion 

We have seen that various moments of branch current give rise to different exponents, namely expo- 
nent b associated with the noise amplitude, /3l associated with resistance scaling and d associated 
with the mass of the fractal. The relation between D(A) obtained in Section 3 and the exponent b 
for normalised noise is as follows: 

|§ ~ PR ~ L-». (19) 



Now, S Rr ~ J2 P \ j p 



This gives: 
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or, 



Now, 



— ^ Pr t ~ L r . (20) 



^ - X(PR r ) ~ 2- b (21) 

PRr 

b = MVA^)) (22) 

In 2 v ' 

KPR) = (23) 



+ 1 



or, 



A r \\{R) 

Substituting this in cquation(8) gives: 



" = Si (26 > 

This expression gives the respective values of the exponent b as 1.1844, 1.0629 and 0.9269 for 
the 3, 4 and 5-simplex respectively. It is clear that the inequality d > b > — /3l is satisfied in each 
of the three cases. 

In the limit n goes to infinity A(2) = (n + 2)/n goes to 1. This is due to the fact that a large 
number of parallel equi-resistance paths are available for current flow. Such a large number of paths 
are available that in going from one order to the next we are in effect not altering the equivalent 
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resistance. The exponent @l decreases in magnitude as we go to higher dimension, implying that 
resistance becomes less dependent on the length of the fractal. 

With regard to flicker noise, we have seen that the scaling relation becomes increasingly complex 
as the order of simplex is increased. The noise versus resistance exponent Q is defined by the 
following: 

Q = 2 +l (27) 

where t and k are given by: 

i2~(Ap)~* (28) 

and 

§ ~ (A P )- fe . (29) 

The experimental measurements [19-20] of t and k were made on 2d-carbon-vax mixtures and 
found to be 2.3 ± 0.4 and 5 ± 1 respectively. The direct plot of Sr versus R leads to S r ~ R® where 
Q = 3.7 ± 0.2. The value we obtain for Q is in agreement with this as is clear from Table 1. 

However, similar measurements on two dimensional films and metallic films have given values of Q 
differing from what we predict. Perhaps instead of taking the n-simplex lattice, if one considered a 2- 
d Sierpinski gasket [13-14] better results could be expected. For all n > a, there is a finite dimension 
matrix whose largest eigenvalue will give the characteristic exponents. The matrix elements are 
function of n and hence eigenvalues will be the well defined function of n. But for n < a the result 
will be obtained by smaller matrix. The generalization to higher value of a and rescaling factor 
b > 2 is under progress. 
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Tabic 1: Exponents calculated for 3-, 4- and 5-simplex lattices. 



simplex 


d 


b 




Q 


3 


1.585 


1.184 


0.737 


3.607 


4 


2.000 


1.063 


0.585 


3.817 


5 


3.322 


0.927 


0.485 


3.909 
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